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Abstract 

> . 

in . 

We present a fast hadronization model for the constituent quark plasma (cqp) 
produced in relativistic heavy ion collisions at SPS. The model is based on rate 
equations and on an equation of state inspired by the string phenomenology. This 
e.o.s. has a confining character. We display the time evolution of the relevant phys- 
ical quantities during the hadronization process and the final hadron multiplicities. 

*^ The results indicate that the hadronization of cqp is fast. 

Oh! 

^ : 1 INTRODUCTION 



X 



The production and investigation of quark matter is the main motivation of the high 
energy heavy ion research. In the early days of the history of this search it was assumed 
that the quark matter formed in high energy heavy ion collisions is similar to a quasi- 
stationary quark gluon plasma, consisting of massless gluons and quarks and antiquarks 
with current mass [p], 0, [J. However, as more and more theoretical and experimental re- 
sults were gathered, gradually the picture of a quark matter emerged, containing effective 
propagators and interaction vertices ||. In the matter formed in the CERN SPS 

heavy ion experiments the effective mass of the interacting gluons is larger than that of 
the dressed quarks 0. 

On the ground of these theoretical indications we expect that in the CERN heavy ion 
experiments not an ideal quark gluon plasma, but a constituent quark antiquark plasma, 
cqp, is formed in some intermediate state of the reaction 0. || . 

*Talk given at the 4th Int. Conference on Strangeness in Quark Matter, July 20-24, 1998, Padova, 
Italy 
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2 QUARK MATTER WITH STRING-LIKE INTER- 
ACTION 



In our model we assume that at the beginning of the hadronization the matter consists 
of massive quarks and anti-quarks. In the time evolution of the system quarks and anti- 
quarks form diquarks, anti-diquarks, mesons, baryons and anti-baryons. We assume that 
the mixture of all of these particles is in thermal equilibrium which can be characterized 
by a temperature. For the representation of the interaction of the colored particles we 
introduce an extra term into the free energy, which is inspired by the string picture. 

As a starting point we consider a mixture of ideal gases of massive quarks, diquarks, 
mesons and baryons, and their respective anti-particles. The corresponding free energy is 

F id = J2 T ^ ( ln ^- x ) + E™^> (!) 

i \ *>th / i 

with Maxwell-Boltzmann statistics for the non-relativistic massive matter: 

N * h = Vd *Jj^ e ~ p2/2miT - ® 

Here d; L = (2sj + l)q are spin and color degeneracy factors. The chemical potentials in 
an ideal gas mixture are 

/Xi,id = Tm T7 — + m i- (3) 



N i>th 



The total entropy is given by 



;,th / 

The energy and pressure of such an ideal, non-relativistic mixture is given by 



ii 



^ = y:n~. (5) 



The energy conservation for the expanding system is expressed as 



dE + pdV = 0. (6) 

At this point we note an important difference to the canonical approach to color 
confinement transition: in an (ideal) mixture of quarks and hadrons the occupied volume, 
V, is the same for both components, V q = Vh = V, the pressure contributions p q and ph 
are additive. On the other hand in the application of the Gibbs criteria of a phase co- 
existence the volumes V q and 14 are additive V q + Vh = constant, and the partial pressures 
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are equal, p q = Ph = P, in phase-equilibrium. In our physical picture of hadronization 
there is no phase coexistence and the Gibbs criteria do not apply. Colored particles 
and color-neutral clusters, pre-hadrons are distributed in a common reaction volume and 
chemical reactions convert eventually the quark matter into a pure hadronic matter. 



The expansion law follows from eq.(||) and eq. (j^) as being 

£ m t N t + £ N t + J2 + T^Y, *i = 0- (7) 

% i i i 

Due to the foregoing hadronization the number of particles decreases, N{ < 0). there- 
fore this process re-heats the system. Cooling effects are due to the expansion (V/V = <9 M u M > 0) 
and rest mass creation (J2i m iNi > 0). 

In order to ensure the color confinement we supplement the model by the following 
confinement principle: all particles carrying color charge (quarks, diquarks, anti-quarks 
and anti-diquarks) will be penalized by a free energy contribution stemming from strings. 
The number of strings is proportional to a weighted sum of the number of color charges, 

Q = ]T qi Ni, (8) 

i 

Here % = for hadrons, = 1/2 for quarks and anti-quarks, and = 3/4 for diquarks and 
anti-diquarks. The higher effective charge of diquarks reflects a possibly higher number 
of in-medium partners, to which a string is stretched. 

The average length L of a string depends on the density of colored objects, as 
L = n^ 1 / 3 , where 

n c = J2N t /V, (9) 

i=c 

and the summation i = c excludes color neutral particles (hadrons) . So the free energy of 
the ideal quark matter - hadron matter mixture, F = + AF, is supplemented by the 
following contribution of strings: 

AF = a s n^ 3 Q, (10) 

with the effective string tension a s ~ 1.0 GeV/fm. 

This additional free energy comprises the non-ideality of the equation of state we 
use. Since this addition is proportional to the volume V and the rest depends on densities 
only, it satisfies thermodynamical consistency requirements due to its construction. 

While there is no new contribution to the entropy, S = S^, the pressure, the energy 
and the chemical potentials of colored (g^ ^ 0) particles receive important modifications: 



P = Pid - 



1 -1/3 Q 

E = E id + a s n: 1/3 Q, 

Hi = fi iM + a s n~ 1/3 (qi- -q), (11) 
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with q = Q/(Vn c ). Hadronic chemical potentials have no modifications at all. 

This non-ideal completion of the equation of state influences both the expansion and 
cooling and the changes of particle composition. The correction to eq.(f|) is due to the 
non-ideal chemical potentials for colored particles. The first principle of thermodynamics 
for this system has the form 



TdS = TdS id = (dE + pdV - ]T fudNi j , (12) 

V i ) id 

the non-ideal cooling law becomes 

(dE + P dV) id + £(/%- ^-, id ) dNi = . (13) 

i 

Accordingly eq.(|7|) is supplemented by a new generic term due to the non-ideal equation 
of state, 



(14) 



f 2 V 2 S i (m i /T)iV i - fi iM /T)N t , 

T 3V E.iV, 3 3 E.iV, 

The additional term, 

E^iW s £^V3 ( (15) 

i i=c 

is negative if color charges become eliminated from the mixture. Therefore color confine- 
ment, causing an extra suppression of equilibrium numbers of quarks and alike particles, 
re-heats the expanding fireball as well as the "normal" chemistry of the ideal quark - 
hadron mixture. 

The only physical effect besides a fast expansion - which has, however, kinematical 
limits stemming from scaling relativistic expansion - that can cool the mixture sufficiently 
is rest-mass production. 

Since in our equation of state we have an explicit interaction energy between the 
quarks, our effective quark masses should be less than that given in Ref.0. We shall use 

the following values: (m u = [^«o + m t/J 1 ^ 2 > m d = [ m do + m t/J 1//2 > m s — [ m so + m ??t] 1/ ' 2 ; with 
thermal mass m th = 0.15 GeV and m u0 = m d0 rs 0, m s0 = 0.15 GeV. The clusters have 
a mass according to the average mass extra to the summed valence quark masses of the 
two lowest lying hadron multiplets: the pseudo-scalar and vector meson nonets and the 
baryon octet and decuplet, respectively (Am mcs = Am,di q = 0.408 GeV, Amb ar = 0.805 
GeV). 



2.1 Initial state 

The initial energy density — distributed along the beam direction between —To sinhrfa and 
To sinh r)o — can be related to the center of mass bombarding energy y/s in the experiment, 

£ ° = TTi^rfV (16) 
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On the other hand the initial invariant volume dual to dr at constant r = r is given by 



V = nB%To2rio. (17) 

The initial internal energy (i.e. the energy without the collective flow of a fluid cell) at 
r = To is therefore less than a/s for finite 770: 

E = e V = -^-Vi. (18) 
shr/o 

At the CERN SPS experiment 770 ~ 1.85 (due to some stopping), R « 7 fm, r ~ 0.5 
fm and we obtain Vo ~ 431 fm 3 and Eq « 2.13 TeV. Compared to the total energy 
of about yfs = 3.4 TeV (carried by about 390 participating nucleons in a central Pb- 
Pb collision) approximately two third of the energy is invested into rest mass of newly 
produced particles and thermal motion and and one third into the flow. 

Comparing this with an alternative expression for the thermal energy of an ideal 
massive quark matter, 

e = E^(°)K + | t )' ( 19 ) 

one can estimate the initial temperature at the beginning of hadronization. Using our 



standard values for the incoming quark numbers JlO N u (0) = 544, Nd(0) = 626, further, 
assuming that 400 uu, 400 dd and, with f s = 0.22, 176 ss quark anti-quark pairs are 
created in one central collision, we get from the above equation T = 0.18 GeV. We use 
these numbers for the newly produced quark pairs in order to arrive at the experimentally 
measured hadron and strange particle numbers. 



2.2 Hadronization processes 

Chemical equilibrium is not supposed initially, rather a definite over-saturation of quarks 
in the reaction volume. The initially missing color-neutral hadron states - mesons and 
baryons - are formed due to quark fusion processes in a non-relativistic Coulomb potential. 
The rates for different flavor compositions differ mainly due to the different reduced masses 
of quark anti-quark or quark diquark pairs. First of all this influences the Bohr radius 
in the Coulomb potential J7|. Of course, the presence of a medium - which establishes 
the necessary momentum balance after the fusion - also influences the hadronization rate. 
The cross section for such a 2 — > 1 process in medium is 



a 



p\ 3 16M 2 0Fa 



2 



(p2 + l/a?f 



(20) 



with a = l/(am) Bohr radius of the Is state in the Coulomb potential and p is the Debye 
screening length . Here p is the relative momentum of the hadronizing precursors, m is 
their reduced mass and M is the total mass. For a a running coupling constant was used. 

The relative momenta are taken from a random Gaussian distribution, 

dV(p) oc e- p2/2mT d 3 p, (21) 
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at temperature T and reduced mass m. This method allows us to simulate thermally 
averaged hadronization rates, 

R=(<M). (22) 
m 

2.3 Dynamical Effects of Confinement 

Besides the "confining correction" in the equation of state we also apply a dynamic 
confinement mechanism in our model: the medium screening length p occurring in the 
hadronization cross section will be related to strings pulled by color charges trying to 
leave the reaction zone. This way the screening length p is increased as the color density 
decreases: we keep, however, the product p 3 N c /V constant, 

p(t)=p(0)*(N c (0)/N c (t)f 3 . (23) 

Here N c stands for the number of colored objects. 



2.4 Reaction network 

What remains to specify the model is the system of rate equations describing the trans- 
formation of quark matter into hadronic matter. We consider N F = 3 light quark flavors 
u, d and s. 

There are N f (Nf + l)/2 = 6 possible diquark flavors and the same number of anti- 
diquark flavors. The number of quark anti-quark flavor combinations is N F = 9 while 
that of quark or anti-quark triplet combinations is N F (N F + l)(Np + 2)/6 = 10. In the 
hadronizing quark matter we deal with altogether 2*3 + 2*6 + 9 + 2*10 = 47 sorts of 
particles. 

Let us generally denote quarks by Q, diquarks by D, mesons by M and baryons by 
B. The possible fusion reactions are: Q+Q — > D,Q + Q — > D,Q + Q — > M, Q + D — > 
B,Q + U — > B. 

Our model is completed by the system of rate equations. Considering a general 
reaction of type i + j — > k we account for the changes 

dNi = dNj = -Adt, dN k = +Adt, (24) 

cumulatively in each reaction. Here 

A = Rij^NiNj (l - (25) 

with a thermally averaged rate Rij^k (cf. eq.(^)). The changes stemming from different 
reactions accumulate to a total change of each particle sort in a time-step dt. 
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2.5 Chemical equilibrium 



Chemical equilibrium is defined by the requirement that all chemical rates vanish (cf. 
eq. p5|) . It leads to relations like 

/4 q + A*? = ft ( 26 ) 

for each reaction channel. The correspondence between equilibrium chemical potentials 
and equilibrium number densities is, however, in the general case not as simple as for a 
mixture of ideal gases. From eq.flTTD we obtain an implicit equation for the equilibrium 
numbers, 

Nr = N hth e^ t (27) 

with 

(28) 



We call the attention to the fact that this equation does not lead to a chemical equilibrium 
state below a critical temperature. 

Applying eq.(|ll]), the non-equilibrium chemical potentials, and hence the essential 
factors, e _Mi//T , in the chemical rates (eq.(^)) can be expressed as 

e -„/T = e -, ?/ T . e M-2^0 iVf (29) 

In the combinations appearing in the detailed balance factor of the rate equations using 
eq.(|2T)|) we obtain 

1 - e ^ = 1 -lt^ e ^^' <30) 

with 

A// i = 6 i (iV)-6 i (iV^). (31) 

The corrections bi(N) in the non-equilibrium chemical potentials may in general depend 
on the number densities of several other components on the mixture. 

At this point we note that the extra e~ Afl ^ T factors occur for any non-ideal equation 
of state where the correction to the free energy density is a nonlinear function of the 
number densities. 



2.6 Hadronic decays 

The set of rate equations describes the time evolution of the number of all involved 
quark and antiquark clusters. In order to get the final hadron numbers we integrate 
these equations until the number of colored particles becomes negligible. At this time 
the number of clusters are divided between the corresponding hadrons according to the 
spin degeneracies of the multiplets. This way we obtain a number of hadronic resonances 
(in the present version the vector meson octet and baryon decuplet). Finally hadronic 
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decays are taken into account with the dominant branching ratios obtained from Particle 
Data Table JTTJ] . We assume that secondary hadron-hadron interactions have a negligible 
effect on the finally observed hadronic composition. The time evolution of the entropy 
and temperature is obtained by simultaneous integration of eqs.@ and ([7]). 



Numerical results and discussion 



For the parameters describing the initial state we used those given in subsection [2.1|, while 
for the parameters determining the dynamics of the hadronization we used the following 
values: p = 0.2/m, a = 1.4. 




Figure 1. The time evolution of the temperature T (a), the entropy S 
(b), the pressure p (c) and the energy density e (d) of the system together 
with the partial contributions of quarks and hadrons. 

From Fig. la we can see, that at the beginning of the hadronization there is a rapid 
decrease in the temperature due to the rest mass formation of the hadrons. Shortly after 
that, the re-heating starts as an effect of color confinement (see eqs.(28) and (32)). Fig. lb 
shows, that the total entropy is monotonically increasing during the hadronization. In 
Fig.lc one can observe an interesting pattern in the time evolution of the pressure. The 
partial pressure of the interacting cqp rapidly decreases as the number of quarks decrease. 
As the color density drops, this pressure becomes even negative. The increasing hadron 
partial pressure, however, over-compensates this negative value. The partial and total 
internal energy evolution, displayed in Fig. Id, shows, that the hadronization is completed 
at 2 fm/c after the beginning of the process. The decrease of the internal energy is 
compensated by the work of pressure while making the flow. 
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Fig. 2 shows the time evolution of different colored particles and color neutral clusters 
corresponding to different hadrons. The diquarks are produced from the quarks, and then 
they contribute rapidly to the formation of baryons. 




t (fm/c) t (fm/c) 

Figure 2. The time evolution of the number of colored particles and 
color neutral clusters. The line styles of different flavor compositions are 
indicated in the respective figures. 

In the ALCOR model the ratio of hadronic species are determined by the ratio of 
steepnesses of these curves. Since these curves do not cross each other, one can under- 
stand, why the algebraic ALCOR approach to the solution of rate equations is a good 
approximation. 

In Table 1 the hadron numbers obtained with the Transchemistry model and those 
obtained with the ALCOR |10| and the RQMD |L2| models are shown together with 
the few published experimental data. Table 2 shows a comparison for the multi-strange 
baryon ratios. While in many cases there are intriguing agreements, in some other cases 
there are some discrepancies. The reason for this may originate from two sources: i) the 
experimental data referred here are the production ratios in the overlap window of the 
detector acceptances. Thus, if the momentum distribution of the two particle sorts is not 
the same, then these ratios are not equal to the total number ratios, ii) these calculated 
values are more sensitive to the simplifying assumption, that the hadronization happens 
into the lowest energy baryon octet and decuplet and the two lowest energy meson octet. 
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Table 1: Total hadron multiplicities for Pb + Pb collision at bombarding energy 158 
GeV/nucleon. The displayed experimental results are from the NA49 Collaboration: a is 
from ||13|| ; b is from [Q; c is from |j[5| ; * is estimated from {-ft' - } and {-Kl?}. Theoretical 
results are from the Transchemistry, ALCOR and RQMD ("ropes + no re-scattering" 
version) model. Here it is A - like = A + S° + + S° + 
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Table 2. Strange baryon and anti-baryon ratios measured by WA97 Collaboration [16 1 
and obtained from Transchemistry and ALCOR model for Pb + Pb collision at 158 
GeV/nucl bombarding energy. The experimental data are the production ratios in the 
overlap window of the detector acceptance. 



4 Conclusion 



In this paper we presented a new model for the hadronization of a quark antiquark plasma 
(cqp) based on rate equations in a quark matter - hadron matter mixture. The color 
confinement was taken into account by using consistently a plausible equation of state 
motivated by the string model. Our results presented in the Figures clearly show a very 
fast hadronization. Observing the shape of the time evolution of different hadron multi- 
plicities, it became understandable, why the simpler algebraic approximation, applied in 
the ALCOR model, works so well. The comparison with the existing experimental data 
indicate, that it is possible, that in the PbPb collision at SPS a piece of matter is formed, 
inside which the massive quarks and anti-quarks interact with a string like mean field. 

Finally we emphasize, that this type of phenomenological investigations are neces- 
sary, as long as the hadronization of the quark matter as a non-equilibrium, non-static, 
non-perturbative process, cannot be described with other methods. 
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